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Abstract 

We consider the Einstein-Gauss-Bonnet equations in five dimensions including a negative cosmological 
j /^ ■ constant and a Maxwell field. Using an appropriate Ansatz for the metric and for the electromagnetic 

fields, we construct numerically black holes with two equal angular momenta in the two orthogonal 
space-like planes of space-time. Families of such solutions, labeled by the angular momentum and by the 
electric charge are obtained for many representative intervals of the Gauss-Bonnet coupling constant a. 
It is argued that, for fixed values of a, the solutions terminate into extremal black holes at (a-dependent) 
critical values of the angular momentum and/or of the electric charge. The influence of the Gauss-Bonnet 
coupling constant, of the charge and of the cosmological constant on the thermodynamics of the black 
Ch ' holes and on their domain of existence is analyzed. 
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^ ! 1 Introduction 

> 

Q\ ' There has been a huge activity in the study of fundamental interactions in space-times involving more than 

four dimensions in the last years. The standard Einstein-Hilbert action for <i-dimensional space-times (with 
d > 4) is not the only choice leading to second order equations of motion. Suitable combinations of terms 
with higher powers in the curvature can be added such that second order field equations result. In the case 
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d = 5, the most general theory of gravity possessing this feature is Einstein-Gauss-Bonnet (EGB) gravity. 
Higher derivative curvature terms occur in many contexts such as in semi-classical quantum gravity and in 
the effective low-energy action of superstring theories. The Gauss-Bonnet term appears as the first curvature 
stringy correction to general relativity [TJ [5] when assuming that the tension of a string is large as compared 
to the energy scale of other variables. Since black holes seem to be the most fundamental solutions of the 
conventional four-dimensional theory of General Relativity, the construction of such solutions in theories 
formulated in d-dimensions and involving higher derivative curvature terms is a natural generalisation. The 
inclusion of the Gauss-Bonnet term in the theory leads to new features of the solutions (see e.g. [3J 0] for 
reviews and further references). Due to the non-linearity of the equations, however, it is very difficult to find 
non-trivial exact analytical solutions of the field equations with higher derivative terms. The main solution 
existing in a closed form is the Gauss-Bonnet counterpart of the Schwarzschild-Tangherlini black hole [5j [6] . 
Static Anti-de Sitter (AdS) black hole solutions in EGB gravity are also known in closed form presenting a 
number of interesting features (see e.g. [7J [HI [HI 02] and references therein). It is of interest to generalize these 
solutions to the rotating case. In five dimensions, rotating solutions with a negative cosmological constant 
have been considered in [TT] within a perturbative approach. The numerical construction of rotating black 
holes (in 5-dimensional EGB gravity) with equal angular momentum in two orthogonal planes was addressed 
in [T2J E3] an d [E], respectively for asymptotically Anti-de Sitter (AdS), de Sitter (dS) and flat space-time. 
In [121 H3] the emphasis was put on the construction of the action and the counterterms necessary to make 
the variational principle and asymptotic charges well defined. Families of black holes were constructed and 
the physical parameters characterizing them were estimated numerically. A systematic study of the domain 
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of existence of the asymptotically flat black holes with respect to their rotational parameter and the Gauss- 
Bonnet coupling constant is reported in |14j . It turns out that for fixed value of the Gauss-Bonnet coupling 
constant, the rotating solutions exist up to a maximal value of the angular velocity at the horizon. The 
construction of the domain of existence leads to numerical difficulties at the approach of the limiting extremal 
configuration. 

In this paper, we extend the study of the domain of existence of the rotating black holes in the presence 
of a negative cosmological constant (AdS case). For a vanishing cosmological constant, the solutions exist 
for arbitrary values of the Gauss-Bonnet coupling, which we call a in the following. When a negative 
cosmological constant is supplemented to EGB gravity, the solutions exist only on a limited interval of a. 
This occurrence of a maximal value for a is examined within our approach of the field equations. In addition 
to the purely gravitating EGB theory, we further supplement the EGB action by a Maxwell term. Solving 
the corresponding field equations we are able to construct families of charged, rotating black holes of the 
Einstein-Gauss-Bonnet-Maxwell equations. The influence of the electromagnetic charge on the domain of 
existence of the black holes is emphasized. 

In section 2 we give the model, we present the Ansatz and also discuss the boundary conditions. Several 
families of exact solutions are reviewed in section 3. Our results are described in section 4 and illustrated 
by several figures. Section 5 gives a summary and conclusions. 

2 The Model 

2.1 The action and field equations 

We consider the Einstein-Gauss-Bonnet-Maxwell action with a cosmological constant A given by 

where R is the Ricci scalar, F 2 is the standard Maxwell density and Lqb represents the Gauss-Bonnet 
density : 

Lqb = R 2 — 4R/j, v R fll/ + R IJ , V(Tr R ,M " Tr . (2-2) 

In (|2.I[) . G and a denote the Newton constant and the Gauss-Bonnet coupling constant, respectively. It is 
usual to relate the cosmological constant to the AdS radius £ according to A = — (d — 2)(d — l)/(2£ 2 ). The 
notations and conventions are the same as in |12j . 

The variation of the action (|2.1j) with respect to the metric tensor and the electromagnetic fields leads 
to the field equations 

la , 

R M u- -Rg fW + Ag fW + -H^=T^ , d^V^gF^) = (2.3) 

where is the electromagnetic stress-momentum tensor while the Gauss-Bonnet tensor is given by 

Hfj,u — 2(R l _ iaKT R l/ ah ' T — 2R iipva R pa — 2R lLa R a v + RR^) — -Lgb9^v ■ (2.4) 

For a well-defined variational principle, one has to supplement the action (|2.ip with the Gibbons-Hawking 
surface term 

IP = ~ I d^x^K , (2.5) 
onG J sm 

and its counterpart for Gauss-Bonnet gravity [5] 

lf B) =-7Trf (J - 2 Ga^) , (2.6) 
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where 7^ is the induced metric on the boundary, K is the trace of the extrinsic curvature of the boundary, 
G a b is the Einstein tensor of the metric 7 (, and J is the trace of the tensor 

J ab = X -(2KK ac K c b + K cd K cd K ab - 2K ac K cd K db - K 2 K ab ) . (2.7) 

2.2 The Ansatz 

While the most general EGB rotating black holes would possess two independent angular momenta and a 
more general topology of the event horizon, we restrict ourselves here to configurations with equal magnitude 
angular momenta and a spherical horizon topology. The suitable metric Ansatz reads |16j 

ds2 = -7TT + g(r)d9 2 + h(r) sin 2 9 (dtp - w(r)dtf + Mr) cos 2 9 (dip - w(r)dtf (2.8) 

fyr) 

+ (g(r) - h(r)) sin 2 6 cos 2 9 {dtp - #) 2 - b{r)dt 2 , 

where 9 £ [0,7r/2], (tp,ip) € [0, 2ir], r and t being the radial and time coordinates, respectively. For such 
solutions, the isometry group is enhanced from Rx U(l) 2 to Rx U{2), where R refers to the time translations. 
This symmetry enhancement allows to transform the Einstein-Gauss-Bonnet equations into a system of 
ordinary differential equations (ODEs). In our construction of the solutions, we fix the residual freedom 
on the coordinate r by taking g(r) = r 2 . The variable r defined in this way is the counterpart of the 
Schwarzschild radial coordinate. Note that in several papers dealing with black holes in higher dimensions 
(see e.g. [TB]) isotropic coordinates are adopted. 

When the Maxwell term is involved, the ansatz (|2.8I) has to be supplemented with a consistent 
parametrization of the electromagnetic field. The vector field is parametrized in terms of two functions 
according to 

A^dx" = V{r)dt + A{r) (sin 2 9dp + cos 2 9di})) (2.9) 

where V{r), A{r) represent respectively the electric and magnetic potentials. The ansatz (|2.8|) . (I2.9[) trans- 
forms the set of Einstein-Gauss-Bonnet-Maxwell equations into a system of six differential equations for the 
radial functions f{r),b{r),h{r),w{r) and V(r), A{r). These equations can be combined in such a way that 
the equation corresponding to the metric function f{r) is of first order. The equation corresponding to the 
electric potential depends only on dV/dr and can be integrated in terms of the other functions. For all the 
other functions the equations are of second order. The form of these equations is involved and not given 
here; it can be found e.g. in [TUfTS]. In the following the prime will denote the derivative with respect to r. 

2.3 Boundary conditions 

We first discuss the regularity properties of the metric fields at the horizon (the conditions for the electro- 
magnetic potentials will be given afterwards). We are interested in solutions possessing a regular horizon at 
r = Th- This requires the conditions f{rh) — 0, b{rh) = 0. The fields can be expanded around the horizon 
leading to 

f{r) = fx{r-r H ) + 0{r-r h f , b(r) = h{r - r h ) + 0{r - r h ) 2 , (2.10) 

h{r) = h +O{r-r h ), w(r) = w h + Wl {r - r h ) + 0{r - r h ) 2 , (2.11) 

where f±, b\, ho are constants which are determined numerically while Wh = w{rh) (we also use the notation 
17 = Wh) represents the angular velocity at the horizon. This parameter controls the rotation of the black 
hole and has to be given by hand. The regularity of the fields at the horizon leads to an extra condition, say 
Ti{f,b,b',h,h',w,w') r — rh = 0, which needs to be imposed. The polynomial Ti is lengthy and is not given 
explicitly here. 
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In the asymptotic region, the metric is assumed to be asymptotically flat if A = or Anti-de Sitter if 
A < 0. The dominant terms of the asymptotic expansions of the fields are (see [T2"] ) 

f(r) = £ (l - f^pj + 1 + + 0(l/r 4 ) , 6(r) = £ (l - f^j + 1 + ^ + 0(1/0 (2.12) 

fc( r) = r 2 + ^ 2 ^^ + /l^W) + Wr 8 ) , »(r) = ^ + 0(l/r 8 ) (2.13) 

Imposing the black hole condition, the regularity at the horizon, the asymptotic form of the metric and 
fixing the angular velocity at the horizon through fl leads to a consistent set of seven boundary conditions 
specifying, in principle, a solution once and a are chosen. The formula (|2.12[) reveals that the expansion 
holds for a < I 2 /2. In the limiting case a = £ 2 /2 the asymptotic expansions are different; e.g. the dominant 
term of the rotating function w(r) is 0(1/ r 2 ) and for the metric function h we have h(r) — r 2 +Ho + 0(l/r 2 ). 
This contrasts with (|2.12l) . The corresponding model, called Chern-Simons gravity, was studied in [Tj5] and 
more recently in [2"U] . 

In the presence of an electromagnetic field, the boundary conditions for the functions V(r) and A(r) can 
be fixed as follows : at r = we set V(r^) = and the condition for the fields to be regular at the horizon 
leads to an equation T2(f,b,b',h,h',w,w',V',A,A') r — rh = 0. The polynomial T 2 (which comes out to be 
the same for the equations for V(r) and A(r) and is again too involved to be given explicitly) is independent 
from Ti mentioned above. In the asymptotic region, the electromagnetic potential has the following form 

V(r) = V mf + ^L + 0(iy) , A(r) = ^ + 0(iy) , (2.14) 

where Vi n f, q and q m are constants. In the absence of rotation (i.e. if w(r) = 0) the field equation 
corresponding to the magnetic potential is trivially fulfilled by A(r) — 0. Practically, we impose V'(r 
oo) = —q/r 3 with a given value of the electric charge q. The last boundary condition consists in imposing 
the product r 2 A(r) to approach a constant. The value of this constant, interpreted as the magnetic charge 
(say q m ) of the solution, depends on a, q, and can be extracted from the numerical data. Static solutions 
have Q = q m = 0. 



2.4 Physical quantities 

Physical quantities characterizing the black holes can then be extracted from the knowledge of the different 
fields. In particular, the Hawking temperature Th , the area Ah of the horizon are given by 
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The energy E and the angular momentum J are given by 



Th = ^^ , A H = V 3 rlVh(r^ . (2.15) 



„ V 3 f 2 - Ab 2 2a V 3 2a 

E = 8^G^2— V 1 " ^ ' J = 8^G W4 V 1 " T ■ (2 ' 16) 

In the figures in the next sections Tjj will be given in units of An and E,J in units of Vz/(&nG). The 
electromagnetic properties of charged, rotating solutions can be characterized, namely by the electric charge 
Q, the magnetic moment fi m and the gyromagnetic ratio g 

= 8irG q ' Mm = 8^G qm ' 9 = 2 ^T ■ (2 - 17) 

Another relevant quantity for the analysis of the thermodynamics is the entropy S of the black hole. It 
can be expressed in terms of the horizon and of the value hh = h(rh) of the metric function h(r) at the 
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horizon (see e.g. [12]) according to 

S=^(rlVW^+^VW^^-^)) • (2-18) 

The study of S as a function of the temperature reveals some properties about the stability through the 
specific heat Cn which is given by : 

C n =T H (dS/dT H )n H . (2.19) 



3 Known solutions 

In this section, we review the explicit solutions of the equations which are available in several specific limits. 
The fact that we recover these solutions with our numerical techniques provides a good crosscheck of the 
accuracy and robustness of our method. 

3.1 Myers-Perry solutions 

We set A = —6/£ 2 where £ is the AdS radius. For pure Einstein gravity, the Myers-Perry (MP) black holes 
[2Tj have the form 

2MS , 2Ma 2 , 2 f it 2Ma 2 \ L r 2 f 2Ma 



with S = 1 — a 2 1 ' i 2 . The parameter a is related to the angular momentum and M is related to the horizon 
Th with f(rh) = 0. Setting for simplicity ry l = 1 the family of MP solutions present the following properties : 

, _ 2(2a 2 (l 4 + 2l 2 + l)-l 4 -2l 2 ) _ -2(2a 2 (l 4 + 2l 2 + l)-l 4 -2l 2 ) 
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Figure 2: Some parameters characterizing the rotating solutions for a = 2.0 as functions of J. 

l + £ 2 ^ 4- P 2 



£ 2 ' 2{£ 2 ~a 2 {l+£ 2 )) ' 

From now on we define Q = Wh = w{rh). In particular, the solution is singular for a 2 = £ 2 /(l + £ 2 ). The 
rotating BH are therefore limited to a 2 < £ 2 /(l + 1 2 ). Let A = a 2 - (£ 2 (£ 2 + 2)/2(£ 2 + l) 2 ). For A < the 
solutions have f'(rh) > and for A > they have f'(rh) < and admit a second horizon for r = > 1. 
Since we want to construct the solutions from the event horizon to infinity, it make sense to consider the 
solutions only for A < and for r G [1, oo]. In the limit A = the black hole is extremal. The corresponding 
values of the parameters a, M, say a ex , M ex , are then fixed in terms of the AdS radius £ : 



.2 _^ 2 + 2) „ 1 



a ex ~ 2 (f2 + iy ' Mex - p ■ ( 3 - 1 ) 
The energy and angular momentum of the extremal black holes can further be determined and read 



V 3 6l 4 + 13l 2 + 8 V 3 y/2(£ 2 + 2)(l+£ 2 ) 

ex ~ SirG 2£± ' ex ~ SttG £ 3 ' 1 ' 



3.2 Anti-de Sitter-Reissner-Nordstom solutions 

Static charged black holes are known in an explicit form. In the limit w(r) = A(r) = 0, the field equations 
are solved by 

r 2 9 M n 2 a 

m . K r) = 1 + F - + 2(d _ 2)( ;_ 3) ,. aJ _„ , m - «W - r> , V(r) . , (3.3) 

where q is a parameter related to the electric charge. Setting for definiteness £ 2 = oo and = 1, we find 
M = (12 + q 2 )/24 and Tjj = (2 — g 2 /6)/47r. The black hole is extremal when the charge become large 
enough: i.e. for q — \/12. 
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Figure 3: Dependence of the angular momentum J on the temperature Th for (asymptotically flat) rotating 
black hole with a = 0,0.35,1,2 (solid lines) and the charged counterparts (with Q = 1, dashed lines) 
(left). Dependence of Th on Wh for rotating black hole with a = 0,0.35,1,2 (solid lines) and the charged 
counterparts (with Q — 1, dashed lines) (right). 



3.3 Static Einstein-Gauss-Bonnet solutions 

Concentrating now on the EGB equations and abandoning rotations (i.e. with w(r) = A(r) = 0) a family of 
exact solutions [SJ [SJ [§] [TO] exists with : 



/ = b = 1 



r\ / , 2M 

- 1- Jl + 2a — 
a V r* 



1 



-)) > 9 = h = 



w = 



(3.4) 



£ 2 12r 6 ' 

Setting q = and I = co (i.e. A = 0) shows that these solutions exist for arbitrary values of a. If A < 0, 
they exist for a < Setting = 1, the charged EGB static black holes have 



M 



12 + 6a + q 2 
24 



12 



6(1 + a) £ 2 (1 + a) 



(3.5) 



Knowing that black hole solutions exist in these two particular cases, a = and 57 = we expect that 
rotating black holes of the Einstein-Gauss-Bonnet equations should exist in some domain of the a — O-plane. 
The purpose of the next section is to determine this domain with some accuracy. 



4 Numerical results 

To our knowledge, no explicit solutions exists in the generic case Wh > 0, a > 0. The system of differential 
equations has to be solved numerically or perturbatively (e.g. using the angular momentum as expansion 
parameter, see |11|). We solved the system numerically by using a collocation method |22j for numerous values 
of a and of Wh- For the main part of the analysis we set = 1 without loss of generality since this corresponds 
simply to fixing the scale of the radial variable r. The study of the quantity 6'(r/j,w/j) = ^\ r = rh . w ( rh )= Wh 
(i.e. the value of the derivative db/dr for the solution with w(r^) = Wh at the horizon) is of some interest 
for the understanding of the results. Referring to the case a = and using an argument of continuity, it 
could be expected that a branch of rotating solutions exists for a > 0, stopping at a maximal value w m {ot) 
such that b'(rh,w m ) — 0. This implies that an extremal black hole with Th = is approached. 
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Figure 4: Left :The parameters w m and J m corresponding to the extremal black holes as functions of a for 
several values of I. Right : The mass of the extremal black holes M m as function of a for several values of £ 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 5: Tjj — J diagram for different values of a and of I. 
4.1 Asymptotically flat solutions. 

The solutions corresponding to A = have been studied in [T3], but we find it useful to summarize the 
main features for completeness. Also, the deviation of the pattern of AdS solutions from the asymptotically 
flat solutions will clearly appear. Fixing a > and increasing Wh gradually, rotating black holes can be 
constructed up until a maximal value oiwh, i.e. Wh = w m (a). Profiles of the solutions corresponding to a = 2 
are presented in Fig. Q]for the non-rotating solution Wh — and Wh ~ 0.68 (for a — 2, we have w m ps 0.72). 
The study of the evolution of the parameters b'(rh), f'(fh) as functions of Wh suggests a natural explanation 
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Figure 6: Evolution of some parameters with a for £ — 2 and Wh = 0.9 



for this result and reveals the pattern of the solutions. For Wh = 0, we have b'(rh = 1) = 2/(1 + a), as can be 
computed from the explicit solution (|3.4p . Then, increasing w h , our numerical results show that b'(r h ), f'{r h ) 
both decrease monotonically and reach the value zero in the limit Wh — >• w m . At the same time, the ratio 
f"{ r h)l f'(fh) becomes very large. This critical phenomenon is illustrated in Fig(5]for a — 2.0, where some 
physical quantities characterizing the solutions have been plotted as function of the angular momentum. 
Reiterating this construction for different values of a confirms that the families of rotating solutions of the 
Einstein-Gauss-Bonnet equations systematically approach an extremal black hole for a maximal value of the 
angular velocity Wh ■ It is then natural to study how the physical properties of the black holes are affected by 
the Gauss-Bonnet parameter a. This is illustrated in Fig. [3] : the left part represents a temperature-angular 
momentum plot for four values of a (solid lines). The right part of Fig. [3] shows the temperature as a 
function of Wh- Note that the lines a = are determined analytically by means of 



The different curves show the non-trivial dependence on the solutions on a. In particular w m and J m are 
not monotonic functions of a. For the slowly rotating black holes the Hawking temperature decreases while 
a increases. This statement clearly does not hold for the solutions with large angular momentum. The 
natural question is then to attempt to determine the domain of existence of rotating EGB black holes in the 
a — ^-parameter space. This is summarized in Fig. IU(left), where the values w m are plotted as a function of 
a. The angular momentum of the extremal black hole is also given. This domain turns out to be the region 
below the solid black line of the figure (representing Wh) or equivalently below the solid red line if we consider 
it from the point of view of the angular momentum J. The corresponding energy of the limiting solution 
is plotted in Fig. |4] (right). The effect of the Gauss-Bonnet interaction is clearly perceptible in these plots. 
The analysis was limited to a € [0,2] since the Gauss-Bonnet term is supposed to come out as a correction 
to the Einstein gravity. However, we found families of solutions presenting the same pattern even for larger 
values of a. 
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Figure 7: The entropy S and horizon 77, as functions of the temperature for A = 0, O = 0.2 and three 
different values of a (left). The mass as function of the temperature for A = 0, ft — 0.2 and three different 
values of a (right). 





Figure 8: The entropy S and horizon ru as functions of the temperature for i = 2, fl = 0.2 and three 
different values of a (left). The mass as function of the temperature for £ = 2, fi = 0.2 and three different 
values of a (right). 



4.2 Asymptotically Anti-De Sitter solutions 

We now discuss the asymptotically AdS solutions (i.e. corresponding to I < oo). The numerical results 
show that they approach continuously the asymptotically flat solutions in the limit £ — > oo. Fig. @] contains 
the data for £ = oo, £ = 4 and 1 — 2. The domain of existence of the solutions varies smoothly with the AdS 
radius £. The influence of the cosmological constant on the temperature and angular momentum is apparent 
in Fig. [5j where the curves are reported for two values of a and for £ = 2,4 and £ = oo. As pointed out 
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Figure 9: Profile of a charged, rotating black hole with a — for Q — 1 and Wh — 0.68 (left), and for a = 2 
for Q = 1 and tift = 0.68 (right), respectively. 



already, the main difference between the flat and AdS black holes resides in the fact that the asymptotically 
AdS black holes exist only for a < £ 2 /2, irrespectively of their angular momentum. Constructing a branch of 
solutions with fixed Wh while increasing a reveals that the asymptotic parameters /2,&2 and W4 all diverge 
in the limit a — > £ 2 /2 while the mass and angular momentum (|2.16p remain finite. This is illustrated in Fig. 
[6] In fact the case a — £ 2 /2 corresponds to Chern-Simons gravity [19]. The solutions in this case will be 
discussed elsewhere [25] . 

4.3 Thermodynamics 

Above we have discussed families of solutions with fixed event horizon and examined the influence of 
the other parameters on it. In order to analyze the thermodynamics of the black holes, one usually fixes 
a and £ and constructs the families of solutions obtained by varying r^. Here we emphasize the influence 
of the Gauss-Bonnet constant a and of the AdS radius £ on the pattern of the solutions. In this numerical 
construction, we keep the horizon angular velocity il fixed (more precisely we set £1 = 0.2). The case A = 
was presented in detail in [14]. In Fig. [7] we give the entropy (left) and the mass (right) of asymptotically 
flat black holes corresponding to three different values of a. As pointed out in [14], the non- vanishing Gauss- 
Bonnet coupling constant leads to a branch of stable black holes with Cq > 0. These possess small horizons. 
Note that it turns out to be difficult to solve the equations for — > 0. We hence limit our analysis to 
rh > 0.1. Our results confirm that no stable rotating black holes exist for a = and that increasing a leads 
to the emergence of a stable branch corresponding to small horizons. The branch of stable black hole gets 
larger while a is increased. These results contrast drastically with the corresponding asymptotically AdS 
black holes. The entropy and horizon radius are plotted in Fig. [8]for £ — 2. The presence of a non-vanishing 
cosmological constant reveal a quite different pattern. For a — Fig. [5] reveals that black holes with a 
small horizon are unstable. Increasing r^, there is a phase transition and the solutions corresponding to 
rh > 1.8 are thcrmodynamically stable. Increasing gradually the parameter a reveals that the unstable 
branch corresponding to small black holes merges into a branch of stable solution. Up to three solutions can 
be constructed which have the same angular velocity and same Hawking temperature but different entropy 
and energy. In particular we find two disjoint intervals of possible event horizon values with Cq > 0. The first 
branch occurs for small horizon. In this case we have Cq < (i.e. unstable solutions) for intermediate values 
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of r/i and the specific heat becomes again positive for large . Finally when a gets large enough (typically 
a = 1), the unstable branch disappears and we find only one branch of stable solution. Interestingly, a 
similar thermodynamical pattern was observed in |17j for charged, rotating solutions. Increasing the angular 
momentum in that case leads to the different phases. By contrast, our solutions are not charged but the 
different phases occur by changing the Gauss-Bonnet parameter. The values of the energy of these families 




Figure 10: Estimation of the values ft, J and E corresponding to the extremal black hole as functions of q 
for a = (left). Dependence of the gyromagnetic factor g on SI for several values of Q and of a (right). 

of black holes are given in Fig. [8] (right). In the case when three solutions exist, it shows that the solutions 
with the smallest energy correspond to the smallest horizon. 

4.4 Charged solutions 

Because of the many parameters involved, we limit the construction of charged solutions to asymptotically 
flat solutions (i.e. A = 0) and fix the scale such that = 1. Charged rotating black holes in Einstein-Gauss- 
Bonnet-Maxwell theory can then be constructed for different values of the Gauss-Bonnet coupling constant 
a, of the electric charge q and of the angular velocity at the horizon SI. The explicit solutions available in 
some special limits and the numerical results suggest these solutions to exist only in a limited domain of the 
a, q, parameter-space and to approach extremal black holes on a surface 

S = S(a,q,Q) = . (4.1) 

The determination of this surface represents a formidable task but several limiting curves are available from 
the special cases : 

• The set of points (a, \/l2, 0) belong to E. 

• The critical line (0,q,u max ) is determined numerically and sketched in Fig. [TU] (left). Note that 
numerous solutions of the Einstein-Maxwell equations have been produced, e.g. in [16] . using isotropic 
coordinates. The pattern of solutions might look different from ours since we use a Schwarzschild-like 
radial coordinate and fix r/,. = 1. 

• The critical line (a, 0,uj max ) was obtained first in [14] and reproduced in Fig|4]given by the black-solid 
line on the left. 
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The surface E can be expected to smoothly extrapolate between these three limiting critical curves. Once 
the suitable boundary conditions are implemented, the equations can be solved for generic values of a,tt,q. 
For non-rotating black holes, imposing a charge by means of q =/= still leads to a zero magnetic potential 
A(r) = 0. Charged-rotating black holes necessarily have A(r) 7^ and are further characterized by a 
magnetic charge. We produced indeed several families of charged, rotating black holes. Fixing the electric 
charge q, our results show that the branch of rotating charged black holes exists and possesses the same 
qualitative pattern as that of the uncharged solutions. Profiles of charged-rotating black holes for a = and 
a = 2 are presented in Fig. [5J Physical data characterizing families of rotating, charged solutions charge for 
several values of a are given in Fig. [3] The solid lines corresponding to q = have been mentioned above. 
The effect of the charge (we set q = 1) is presented by the dashed lines for the cases a = and a — 2 only. 
This figure reveals that the presence of the charge changes the domain of existence only little as compared 
to the uncharged black holes. The numerical construction of several other branches (e.g. with a, ft fixed 
and varying q) confirms this tendency and suggests the surface £ to be smooth. Finally, we studied some 
electromagnetic properties of the solutions. In particular, we present the response of the magnetic moment 
and gyromagnetic factor g (see (|2.17|) ) on a change of the electric charge and of the Gauss-Bonnet coupling 
constant. This is summarized in Fig. [TU] (right). It shows that the Gauss-Bonnet term has a tendency to 
increase the factor g. As expected [TS], the gyromagnetic factor tends to g — d — 2 (i.e. g = 3) in the static 
limit. Note that the different curves on this figure stop at the approach of the extremal black hole. 

5 Conclusions 

In this note we have presented arguments supporting the fact that equal angular momentum rotating black 
holes of the Einstein-Gauss-Bonnet equations in five dimensions exist for generic values of the Gauss-Bonnet 
coupling constant. In the case of a negative cosmological constant, the asymptotically AdS solutions exist 
up to a maximal value of the Gauss-Bonnet coupling parameter, i.e. for a < a cr — £ 2 /3. For a — a cr the 
model corresponds to Chern-Simons gravity [19] - The solutions obey a different asymptotic expansion than 
for generic values of a. We were able to trace back this phenomenon with our numerical construction and 
obtain robust solutions even for a — a cr by implementing the appropriate boundary conditions. 
The effect of the Gauss-Bonnet term would be worth studying for other kinds of solutions of the underlying 
Einstein equations with a negative cosmological constant. It would be challenging, in particular, to study 
the response of the NUT-charged solutions constructed in [24] to a Gauss-Bonnet term. Further applications 
of the Gauss-Bonnet black holes to the context of brane gravity could be looked for by following, e.g. the 
ideas of [25] . 

Supplementing the EGB Lagrangian with a Maxwell term leads naturally to charged black holes. The so- 
lutions terminate into extremal black hole at some maximal value of the angular momentum and/or of the 
electric charge. We believe that the numerical estimates of the domain of existence of charged-rotating black 
holes in EGB gravity (which we studied in the case A = only) is correct qualitatively and quantitatively. 
This could be confirmed independently by using an appropriate attractor technique along the lines of [26] . 
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